1 



The Theory of Nuclear Forces: 
^ ; Is the Never-Ending Story Coming to an End? 
O 

(N ! R. Machleidt 
c3 ' a Department of Physics, University of Idaho, Moscow, Idaho 83844, USA 

; 

I review the recent progress in our understanding of nuclear forces in terms of an effective 
field theory (EFT) for low-energy QCD and put this progress into historical perspective. 
This is followed by an assessment of the current status of EFT based nuclear potentials. 
^ . In concluding, I will summarize some unresolved issues. 

1. Historical Perspective 

O 



X 



The theory of nuclear forces has a long history. Based upon the Yukawa idea [ [T] and 
the discovery of the pion, the 1950's became the first period of "pion theories". These, 
however, resulted in failure — for reasons we understand today: pion dynamics is ruled by 
chiral symmetry, a constraint that was not realized in the theories of the 1950's. 

The 1960's and 70's represent the main period for theories that also include heavy 
mesons ("meson theories") [E], but the work on meson models continued all the way into 
the 1990 's when the family of the so-called high-precision NN potentials was developed. 
This family includes the Nijm-I, Nijm-II, and Reid93 potentials [[3], the Argonne Vis [H], 
and the CD-Bonn potential [[5j[6]. Later, also the highly non-local potential by Doleschall 
et al. [ 7\ and the "CD-Bonn + A" model by Deltuva et al. [ 8J joined the club. All these 
potentials have in common that they are charge-dependent, use about 40-50 parameters, 
and reproduce the 1992 Nijmegen NN data base with a x 2 /datum » 1. 

Over the past ten years, the high-precision potentials have been applied intensively in 
exact few-body calculations and miscroscopic nuclear structure theory. Already the first 
few applications of the high-precision potentials in three-nucleon reactions [ |9] clearly 
revealed sizable differences between the predictions from NN potentials with a x 2 /datum 
~ 1 (i.e., the new generation of potentials) and NN potentials with a x 2 /datum « 2 (the 
old generation of the 1970's and 80's which includes the old Nijmegen, the Paris, and the 
old Bonn potentials). Thus, once for all, the standard of precision was established that 
must be met by any future work in microscopic nuclear structure and exact few-body 
calculations: The input NN potential must reproduce the NN data with a x 2 /datum « 1 
or the uncertainty in the predictions will make it impossible to draw reliable conclusions. 

In spite of these great practical achievements, the high-precision potentials cannot be 
the end of the story, because they are all phenomenological in nature. Ultimately, we 
need potentials that are based on proper theory and yield quantitative results. 



*Plcnary talk presented at the 18th International Conference on Few-Body Problems in Physics, August 
21-25, 2006, Santos, SP, Brazil; to be published in Nucl. Phys. A. 
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Table 1 

X 2 /datum for the reproduction of the 1999 np database by families of np potentials at NLO and 
NNLO constructed by the Bochum/Juelich group [115]. 





# of np 


Bochum/Juelich 


Bin (MeV) 


data 


NLO 


NNLO 


0-100 


1058 


4-5 


1.4-1.9 


100-190 


501 


77-121 


12-32 


190-290 


843 


140-220 


25-69 


0-290 


2402 


67-105 


12-27 



Since the nuclear force is a manifestation of strong interactions, any serious derivation 
has to start from quantum chromodynamics (QCD). However, the well-known problem 
with QCD is that it is non-perturbative in the low-energy regime characteristic for nuclear 
physics. For many years, this fact was perceived as a major obstacle for a derivation of the 
nuclear force from QCD — impossible to overcome except by lattice QCD. The effective 
field theory (EFT) concept has shown the way out of this dilemma. One has to realize 
that the scenario of low-energy QCD is characterized by pions and nucleons interacting 
via a force governed by spontaneously broken approximate chiral symmetry. 

Based upon this EFT, a systematic expansion in terms of (Q/A x ) u can be developed, 
where Q denotes a momentum or pion mass, A x fa 1 GeV is the chiral symmetry breaking 
scale, and v > [ [TO] . This has become known as chiral perturbation theory (ChPT). In 
contrast to the pion theories of the 1950 's, this new chirally constrained "pion theory" is 
quite successful. Starting with the pioneering work by Ordonez, Ray, and van Kolck in 
the early 1990's [ [TT]. NN potentials have been constructed in the framework of ChPT. 
In the next section, I will discuss the latest status of this work. 

A remarkable fact of ChPT is that it makes specific predictions also for many-body 
forces. For a given order of ChPT, two-nucleon forces (2NF), three-nucleon forces (3NF), 
. . . are generated on the same footing (see Ref. [H2] for a concise review). At leading order 
(LO), there are no 3NF, and at next-to-leading order (NLO), all 3NF terms cancel [ \10\ 
IT3] . However, at next-to- next-to-leading order (NNLO) and higher orders, well-defined, 
nonvanishing 3NF occur [[131 EH]- Since 3NF show up for the first time at NNLO, they 
are weak. Four-nucleon forces (4NF) occur first at N 3 LO and are, therefore, even weaker. 

2. Current status of chiral NN potentials 

As pointed out in the previous section, the research conducted during the high-precision 
period of the 1990's [ H] systematically tested and firmly established the standards of 
precision that must be met by NN potentials to yield sufficiently accurate predictions 
when applied in microscopic calculations of few-body systems. Since in nuclear EFT we 
are dealing with a perturbative expansion (ChPT), the question is to what order of ChPT 
we have to go to obtain the precision that meets those standards. To discuss this issue 
on firm grounds, I show in Table 1 the \ 2 1 'datum for the fit of the world np data below 
290 MeV for a family of np potentials at NLO and NNLO. The NLO potentials produce 
the horrendous xV^atum between 67 and 105, and the NNLO are between 12 and 27. 
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Table 2 

X 2 /datum for the reproduction of the 1999 np database by various np potentials. Num- 
bers in parentheses denote cutoff parameters in units of MeV. 







Idaho 


Bochum/ Juelich 


Argonne 


Bin (MeV) 


# of np 


N 3 LO [EH] 


N 3 LO [Q35] 


^18 [HI 




data 


(500-600) 


(600/700-450/500) 




0-100 


1058 


1.0-1.1 


1.0-1.1 


0.95 


100-190 


501 


1.1-1.2 


1.3-1.8 


1.10 


190-290 


843 


1.2-1.4 


2.8-20.0 


1.11 


0-290 


2402 


1.1-1.3 


1.7-7.9 


1.04 



The rate of improvement from one order to the other is very encouraging, but the quality 
of the reproduction of the np data at NLO and NNLO is obviously totally insufficient 
for reliable predictions. In the literature, one can find calculations employing NLO and 
NNLO potentials in few-nucleon systems. Because of the extremely poor accuracy of such 
potentials, no reliable conclusions can be drawn from such calculations, rendering these 
applications essentially useless. 

Based upon these facts, it has been pointed out in 2002 by Entem and Machleidt [ 
[161EL7] that NNLO is insufficient and one has to proceed to N 3 LO. Consequently, the first 
N 3 LO potential was created in 2003 [EH], which showed that at this order a x 2 /datum 
comparable to the high-precision Argonne Vx$ potential can, indeed, be achieved, see 
Table 2. This "Idaho" N 3 LO potential [[18] produces a x 2 /datum = 1.1 for the world 
np data below 290 MeV which compares well with the x 2 /datum = 1.04 by the Argonne 
potential. In 2005, also the Bochum/ Juelich group produced several N 3 LO NN potentials [ 
[19], the best of which fits the np data with a x 2 /datum =1.7 and the worse with a 
X 2 /datum = 7.9 (see Table 2). While 7.9 is clearly unacceptable for any meaningful 
application, a xVdatum of 1.7 is reasonable, although it does not meet the precision 
standard established in the 1990's. 

I turn now to the pp data. Typically, x 2 for pp data are larger than for np because 
of the higher precision of pp data. Thus, the Argonne produces a xVdatum = 1.4 
for the world pp data below 290 MeV and the best Idaho N 3 LO pp potential obtains 
1.5. The fit by the best Bochum/ Juelich N 3 LO pp potential results in a xVdatum = 2.9 

Table 3 

% 2 /datum for the reproduction of the 1999 pp database by various pp potentials. Num- 
bers in parentheses denote cutoff parameters in units of MeV. 







Idaho 


Bochum/ Juelich 


Argonne 


Bin (MeV) 


# of pp 


N 3 LO [rjs] 


N 3 LO [H9] 


Vis [S 




data 


(500-600) 


(600/700-450/500) 




0-100 


795 


1.0-1.7 


1.0-3.8 


1.0 


100-190 


411 


1.5-1.9 


3.5-11.6 


1.3 


190-290 


851 


1.9-2.7 


4.3-44.4 


1.8 


0-290 


2057 


1.5-2.1 


2.9-22.3 


1.4 
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Figure 1. np phase shifts for energies below 300 MeV and partial waves with total angular 
momentum J < 2. The solid and the dashed curves show the phase shifts predicted by the N 3 LO 
potentials constructed by the Idaho [[18] and the Bochum/Juelich [[19] groups, respectively. Solid 
dots and open circles represent the Nijmegen multienergy np phase shift analysis [[20] and the 
GWU/VPI single-energy np analysis SM99 [I21j. respectively. 



which, again, is not quite consistent with the precision standards of the 1990's. The worst 
Bochum/Juelich N 3 LO pp potential produces a xVdatum of 22.3 and is incompatible 
with realiable predictions. 

Phase shifts of np scattering from the best Idaho (solid line) and Bochum/Juelich 
(dashed line) N 3 LO np potentials are shown in Figure 1. The phase shifts confirm what 
the corresponding % 2 's have already revealed. 



3. The low-energy constants 

The two-pion-exchange contribution to the NN interaction depends on the low-energy 
constants (LECs) of the dimension-two Lagrangian, commonly denoted by Cj, and the 
dimension-three LECs, d\ [ [17] . These parameters have been determined in ttN data 
analyses, see column ttN of Table 4. The most influential LEC is C3 and, so, we will focus 
here just on C3. Not surprisingly, the Idaho and Bochum/Juelich N 3 LO potentials apply 
essentially the same value, namely, C3 = —3.3 ±0.1 GeV -1 , which is on the lower side 
but still roughly within one standard deviation of its ttN determination. There exists 
also a determination of this parameter from the NN data by the Nijmegen group. They 
report c 3 = -5.08 ± 0.28 GeV" 1 in their 1999 pp analysis [ 24] and c 3 = -4.78 ± 0.10 
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Table 4 



Low-ener£ 


sy parameters. 


q are given 


in units of GeV 


1 and di in units 






Idaho [I18| 


Boch. /J. [\19\ 


AT' " 

JNijmegen 




_ AT f pjol p)ql 

VTiV [ \ZZ\ Izo] 


1ST3T C\ 

IN L/U pot. 


t\t3t r\ „~+ 
IN IjU pot. 


9HHQ "D"\S/"A r I9£l 

zUUo tr WA [ Izo] 










(1999 rWA [1241 ) 


Cl 


-0.81 ± 0.15 


-0.81 


-0.81 


-0.76 










(-0.76) 


C2 


98 -1- 9 1 } 
o.Zo zt U.Zo 


Z.oU 


o.Zo 




C3 


-4.69 ± 1.34 


-3.20 


-3.40 


(") 

-4.78 ±0.10 










(-5.08 ±0.28) 


C 4 


3.40 ± 0.04 


5.40 


3.40 


3.96 ±0.22 










(4.70 ± 0.70) 


di + d 2 


3.06 ±0.21 


3.06 


3.06 




d 3 


-3.27 ±0.73 


-3.27 


-3.27 




d 5 


0.45 ± 0.42 


0.45 


0.45 




du — di^ 


-5.65 ±0.41 


-5.65 


-5.65 
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Figure 2. 3 F 4 phase shifts cal- 
culated at NNLO using C3 = 
—4.78 GeV -1 . The solid curve 
is based upon the NNLO ampli- 
tude whithout cutoff, while for the 
dashed curve the cutoff used in the 
Nijmegen analysis [[25] is employed. 
Solid dots and open circles like in 
Fig. 1. 



GeV -1 in their 2003 pp + np analysis [ES]- Knowing the sensitivity of NN phase shifts and 
observables to C3, the difference between the NN potential values of about -3.3 GeV -1 and 
the Nijmegen NN analysis value of -4.8 GeV -1 is huge. Off-hand, it is hard to understand 
how two different types of NN analyses can be in such severe disagreement. 

However, the discrepancy is easily explained if one looks into the details of the Ni- 
jmegen analysis [EE]- The Nijmegen group starts from the NN amplitude at NNLO in 
momentum space, Fourier transforms it into r-space to obtain a local potential, and then 
cuts this potential off at r=1.6 fm (i.e., the potential is set to zero for r < 1.6 fm). The 
impact of this cutoff is demonstrated in Figure 2 for the phase shifts, which is a 
representative case. In that figure, the solid line shows the phase shifts derived from the 
NN amplitude at NNLO using c 3 = —4.78 GeV -1 and applying no cutoffs (i.e., the actual 
model- independent amplitude at NNLO is used). The dashed line is obtained when the 
Nijmegen cutoff is applied and, as before, c 3 = —4.78 GeV -1 is used. The difference 
between the two curves in Fig. 2 is entirely due to the cutoff applied in the Nijmegen 



6 



analysis. It is clearly seen that the influence of this cutoff is dramatic and introduces a 
huge systematic error rendering the Nimegen analysis and their results unreliable [ 26J. 

It must be noted that the Idaho and Bochum/ Juelich N 3 LO potentials also use cutoffs, 
but not as harsh ones as in the Nijmegen analysis. 

The best way to determine the LECs from NN data is to apply no cutoffs to the NN 
amplitudes (at NNLO or N 3 LO) in peripheral partial waves and to keep the amplitudes 
of lower partial waves at their empirical values (as determined in phase shift analysis) . In 
such an analysis the value c 3 = —2.7 GeV -1 is obtained [I2"7j. 

4. Conclusions 

The theory of nuclear forces has made great progress since the turn of the millennium. 
Nucleon-nucleon potentials have been developed that are based on proper theory (EFT 
for low-energy QCD) and are of high-precision, at the same time. Moreover, the theory 
generates two- and many-body forces on an equal footing and provides a theoretical 
explanation for the empirically known fact that 2NF 3> 3NF 3> 4NF .... 

At N 3 LO [ [XTl EEE], the accuracy can be achieved that is necessary and sufficient for 
microscopic nuclear structure. First calculations applying the N 3 LO NN potential [HE] in 
the (no-core) shell model [ ESI [2S1 EDI Ell E2] , the coupled cluster formalism [E31 El E3 
1361 137] , and the unitary- model-operator approach [ EE] have produced promising results. 

The 3NF at NNLO is known [ [T3J, [T3] and has been applied in few-nucleon reactions [ 
[HI ESI HQ] as well as the structure of light nuclei [ [4TJ [42]. However, the famous L A y 
puzzle' of nucleon-deuteron scattering is not resolved by the 3NF at NNLO. Thus, the 
most important outstanding issue is the 3NF at N 3 LO, which is under construction [143]. 

Another open question that needs to be settled is whether Weinberg power counting, 
which is applied in all current NN potentials, is consistent. This controversial issue is 
presently being debated in the literature [ HU HU, [46] . 

It may be too early to claim that the never-ending-story is coming to an end, but the 
story seems to be converging — at the same rate as chiral perturbation theory. 
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